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The non-Newtonian behavior of a monodisperse concentrated dispersion of spherical particles was 
investigated using a direct numerical simulation method, that takes into account hydrodynamic in- 
teractions and thermal fluctuations accurately. Simulations were performed under steady shear flow 
with periodic boundary conditions in the three directions. The apparent shear viscosity of the dis- 
persions was calculated at volume fractions ranging from 0.31 to 0.56. Shear-thinning behavior was 
clearly observed at high volume fractions. The low- and high-limiting viscosities were then estimated 
from the apparent viscosity by fitting these data into a semi-empirical formula. Furthermore, the 
short-time motions were examined for Brownian particles fluctuating in concentrated dispersions, 
for which the fluid inertia plays an important role. The mean square displacement was monitored in 
the vorticity direction at several different Peclet numbers and volume fractions so that the particle 
diffusion coefficient is determined from the long-time behavior of the mean square displacement. 
Finally, the relationship between the non-Newtonian viscosity of the dispersions and the structural 
relaxation of the dispersed Brownian particles is examined. 

PACS numbers: 82.70.-y, 82.20.Wt, 47.50.-d, 83.50.Ax 



I. INTRODUCTION 



the Krieger-Dougherty relationship Q : 



The links between macroscopic rheological properties 
and microstructures in colloidal dispersions have been ex- 
tensively investigated for many systems, including disper- 
sions of sterically-stabilized and chargcd-stabilized parti- 
cles in host fluids [J, Q • 

The behavior of a monodisperse dispersion composed 
of solid particles immersed in a Newtonian host fluid 
strongly depends on the volume fraction of the dispersed 
particles $ and the shear rate 7. When the shear rate is 
zero (7 — > 0), the shear viscosity of the dispersion is re- 
ferred to as the zero-shear viscosity 770. In the dilute limit 
(<f> ^ 1), the zero shear viscosity is well approximated by 
Einstein's formula 



770 = ?7(1 + 2.5$), 



(1) 



where 77 is the shear viscosity of the host fluid. In a 
concentrated dispersion, theoretical difficulties become 
rather severe, since the behavior of the dispersed par- 
ticles is complicated by the interactions between the 
particles and thermal fluctuations. In particular, the 
solvent-mediated many-body hydrodynamic interactions 
(HI) between the particles complicate the dynamical be- 
havior. A number of experiments for concentrated dis- 
persions have been performed to reveal the origin of the 
non-Newtonian behavior of these dispersions, and sev- 
eral semi-empirical formulas for r/o have been proposed 
to characterize the experimental results. For example. 
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where [rj\ is the intrinsic viscosity and $m is the pack- 
ing volume fraction at which the viscosity diverges, is 
often used for fitting the experimental data for uniform 
colloidal spheres suspended in non-aqueous media. 

When a dispersion is subjected to shear (7 7^ 0), the 
flow properties of the dispersion show a variety of non- 
Newtonian behaviors, such as shear-thinning and shear- 
thickening. These non-Newtonian behaviors are associ- 
ated with the changing microstructures of the dispersion. 
Several physical mechanisms for these peculiar behaviors 
have been proposed' for example, shear-induced order- 



disorder transitions 
ters of the particles 
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0] , formations of dynamic clus- 
However, a full understanding 
of the relationships between the rheological properties 
and the microstructure has not yet been obtained, de- 
spite extensive studies. 

Computer simulations are very powerful tools in the 
direct investigation of the dynamics of individual parti- 
cles in concentrated dispersions. The Stokesian dynamics 
(SD) method [l^l has been widely used to measure the 
rheology of dispersions and provides valuable information 
regar ding the non-Newtonian behavior of flowing disper- 
sions [14|. [Tsl. [TgI. [iTt. The SD method, however, is based 
on the Stokes approximation, which assumes that all re- 
laxation times associated with fluid motions are short 
as compared with those of the particle, i.e., t^, <^ r^, 
where t^, = pfO? /rj and tb = 2ppa^/9ri. Here, pf and pp 
are the density of the fluid and the particle, respectively, 
and a is the radius of the particle. With this assumption, 
the HI is treated as the Ronte-Prager-Yamakawa (RPY) 
tensor and the lubrication correction. Furthermore, it 
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is assumed that the relaxation time associated with the 
particle's inertia is zero (rg — > 0). These approximations 
are valid for the motion of the particles at time scales 
much greater than the time scales of the relaxation of 
both the fluid and the particle inertia. Therefore, the 
short-time motion of Brownian particles over the kine- 
matic time scale cannot be described by simulation 
methods based on the Stokesian approximation. On the 
kinematic time scale, the dynamic coupling between the 
fluid motion and the particle motion remains strongly. 

Sheared dispersions have a non-dimensional parame- 
ter that includes the particle's inertia, e.g., the parti- 
cle Reynolds number Rep = ppja^/rj = Qtb/'2t^, where 
tai{— 1/7) represents the characteristic time due to shear. 
For most colloidal dispersions, the particle Reynolds 
number is very small. The SD method has a particle 
Reynolds number of zero and, therefore, cannot be ap- 
plied to problems with finite particle Reynolds numbers. 
Typical examples include the motions of dispersions sub- 
jected to strong shear or composed of large particles; in 
these cases, the particle Reynolds number has a relatively 
large finite value. 

Another problem that cannot be treated by the SD 
method is the short-time motion of the particle, in which 
the fluid inertia becomes significant even if the particle 
Reynolds number is very small. For example, the char- 
acteristic time scale is Tj/ ~ 1 (jls for a neutrally buoyant 
particle of 1 /zm in water, and Rcp ^ 10~* for 7 ~ 0.01. 
The effects of the fluid inertia appear as memory effects 
of the particles. For a complete understanding of these 
flowing dispersions, full time-dependent His are required. 

In recent years, several numerical methods have been 
developed in order to accurately simulate dispersions in 
a variety of situations, including those described above. 
These methods of dispersion modeling are based on the 
same approach, which involves resolving the fluid motion 
simultaneously with the particle motion. We refer to 
this approach as the direct numerical simulation (DNS) 
approach. This approach enables us to accurately treat 
the full time-dependent HI. The numerical methods differ 
mainly in the approach used to resolve the His between 
the fluid and particle motions [H, [H, [13, [111, [12, [13, [13, 

In this work, we apply a direct numerical scheme based 
on the smoothed profile method (SPM) [li,!!^,!!!,!!!,!!! 
to a monodisperse concentrated dispersion of repulsive 
and neutrally buoyant Brownian particles in a shear flow. 
In the SPM, the Navier-Stokes equation for the fluid mo- 
tion is discretized on a regular grid, and the Newtonian 
equations for the particle motion are solved simultane- 
ously with the fluid motion. This developed scheme ac- 
counts for thermal fiuctuations, a shear flow, and memory 
effects. The SPM can reproduce the correct short-time 
behavior of a single Brownian particle in a shear flow on 
the kinematic time scale, and the numerical tests for a 
sing le Brownian particle in a shear flow were reported in 

Although some groups have studied flowing disper- 



sions composed of repulsive spherical particles by using 
a DNS method based on the lattice Boltzmann method 
113, all the simulations have been applied to a dis- 
persion of non-Brownian particles and ignore the thermal 
fluctuations of the particles. Systematic analyses have 
not yet fully carried out for concentrated dispersions of 
Brownian particles in a shear flow at flnite Reynolds num- 
bers. 

In the present study, we examine the non-Newtonian 
rheology of concentrated dispersions in shear flow at fi- 
nite Reynolds and Peclct numbers using a DNS method 
that takes into account hydrodynamic interactions and 
thermal fluctuations accurately. Note that the simulated 
situations are different from those of the prevailing exper- 
iments for colloidal dispersions, since even the lowest par- 
ticle Reynolds number in these simulations is several hun- 
dred times larger than those of the experiments. We flrst 
present the simulation method and the manner in which 
the apparent shear viscosity of the dispersions is calcu- 
lated. Three-dimensional simulations are performed with 
periodic boundary conditions, and the non- Newtonian 
behavior of the shear viscosity is obtained for several vol- 
ume fractions. Both the high and low shear limiting vis- 
cosities obtained from the simulations are then compared 
with the Kricger-Doughty relationship. Moreover, we in- 
vestigate the short-time motions of Brownian particles in 
a sheared concentrated dispersion on the kinematic time 
scales. The mean square displacement (MSD) of Brown- 
ian particles is monitored at several Peclet numbers and 
volume fractions, and the long-time diffusion coefficient 
is determined from the long-time behavior of the MSD. 
We also suggest a simple relationship between the non- 
Newtonian viscosity of the dispersions and the structural 
relaxation time of the dispersed Brownian particles. 



II. SIMULATION METHOD 

A direct numerical scheme that implements both a 
shear fiow and thermal fluctuations is briefly explained. 
A more detailed explanation is given in a previous pub- 
lication [sot . We consider a monodisperse dispersion of 
Np repulsive spherical particles of diameter cr in a New- 
tonian host fluid. The particles interact via a truncated 
Lennard-Jones (LJ) potential: 



ULj{r) = 




where r is the distance between two particles and the 
parameter e characterizes the interaction strength. The 
position of the jth dispersed particle is Ri , the transla- 
tional velocity is 1^, and the rotational velocity is il;. 
The time evolution of the ith particle with mass Mi and 
moment of inertia Ii is governed by Newton's equations 
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of motion: 



G 



AT. 



H 



(4) 
(5) 



where F/^ and AT/^ are the hydrodynamic forces and 
torques exerted by the host fluid on the particle. Ff is a 
repulsive force arising from the potential of Eq. ([3]) , which 
prevents the particles from overlapping. GY and Gf are 
random forces and torques, respectively, due to thermal 
fluctuations. These random fluctuations are assumed to 
be Markovian (white or time delta-correlated) and de- 
termine the particle temperature T. The procedure for 
determining the temperature is described in [H, [2^, [13] ■ 

In the SP method, the velocity and pressure fields, 
v{x,t) and p{x,t), are defined on three-dimensional 
Cartesian grids, which consist of fluid and particle do- 
mains. In order to distinguish the particle and fluid do- 
mains on the grids, a smoothed function (j){x,t), which 
is equal to 1 in the particle domains and in the fluid 
domains, is introduced. These domains are separated by 
a thin interfacial domain of thickness ^. The system size 
is [0,LJ X [~Ly/2,Ly/2] X [0,L.]. 

The time evolution of the velocity field is governed 
by the Navier-Stokes equation with the incompressibil- 
ity condition V • t; = 0: 



Pfidtv + V ■ Vv) = V • <T + Pf(f>fp + pff 



shear 



(6) 



where the stress tensor cr = —pi + rjlVv + (Vv)'^}, and 
f^'^"'^ {x,t) is an external force field that is introduced 
to enforce a simple shear flow on the system. The flow 
is imposed in the x direction and the external force is 
introduced as a constraint force so that the velocity field 
satisfies Vx{y) = —jLy/A at y = —Ly/A and Vx{y) = 
jLy/A at y = Ly/A, where y denotes the distance in the 
velocity gradient direction. A simple shear fiow with a 
shear rate of 7 is then approximately produced over a 
range from y = —Ly/A to y = Ly/A. (j)fp represents 
a body force that ensures the rigidity of the particles 
and the appropriate non-slip boundary conditions at the 
fiuid/particle interface, which is further elaborated upon 
in reference [l^, [13] . 

The unit of length is taken to be the lattice spacing 
A, and the unit of time is tq = p/A-^/?/. Unless oth- 
erwise stated, we set A = 1, tq = 1, = 1, p/ = 1, 
Pp = 1, a = 4, e = 1, and a = 8. Assuming disper- 
sions of neutrally buoyant particles of radius l/xm in wa- 
ter at room temperature, our unit length A and time 
To correspond to 0.25/im and 0.0625^s, respectively The 
second-order Runge-Kutta algorithm is used to integrate 
the Newtonian equations. The Navier-Stokes equation is 
discretizcd with a Fourier spectral scheme in space and 
with a second-order Runge-Kutta scheme in time. The 
discretized time step is h = 0.07145. Although this value 
is choosen from the stability condition of the Navier- 
Stokes equation, it can be used safely for the particle's 
equations motion because h is much smaller than the 
Lennard-Jones time unit tm = {Mia^ /e)^^'^ ~ 131. The 



phase diagram of the present (36:18) LJ system depends 
weakly on the system temperature T as well as the vol- 
ume fraction $ of the particles. To avoid the delicate 
issue of crystallization, all the simulations in the present 
paper were carried out at T and $ with which the system 
is in amorphous states. 

To measure the rheological properties of the particle 
dispersion in shear fiow, we can calculate the apparent 
stress cr'^PP of the dispersions in the following manner. 
The momentum equation for the dispersion is formally 
written as: 



dis 



Ptf"'' 

pt = (1 - <l>)pf + <l>Pp, 



(7) 
(8) 



where cr'**'* denotes the stress tensor of the dispersion in- 
cluding the inertia, the pressure and the viscous terms. 
The full stress tensor s of the flowing dispersion is then 
defined by introducing a convectivc momentum-flux ten- 
sor explicitly as: 



PtVV, 



(9) 



where ptvv represents momentum transport by the bulk 
flow of the dispersion. 

Although s cannot be calculated directly, we can ob- 
tain the apparent stress cr°'PP by using the local stress 
s: 



/ dxs 



dx 



[V • (sx)]"^ -xV ■ s 
^ J dx[—xV ■ s] 



1 

V 



dx 



x(ptr'^'^'^^ --(ptv) 



i J dxxptf'^^'^^ - ^ J dxx^^iptv), 



(10) 

(11) 

(12) 
(13) 
(14) 



with a volume V = L^LyLz- In the derivation of Ea. pi[) . 
we use an second rank identity, s = [V • {sx)]^ — xV ■ s. 
In the steady state. 



(|(P..)),= 0. 



(15) 



where Qt denotes time averaging over the steady state. 
The time-averaged apparent shear stress of the dispersion 
can then be written as: 



t = ^{j dxxptt'^'^-^] ^ 



(16) 



Then one can obtain the apparent shear stress under 
steady shear flow from the external force imposed 
in the Navire-Stokes equation. 
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III. RESULTS AND DISCUSSION 

Simulations were performed in a three-dimensional cu- 
bic box, whose side length is L, with periodic bound- 
ary conditions. Most present simulations were done with 
L = 64. where the number of particles are Np — 300, 
400, 450, 500, and 550 for $ = 0.31, 0.41, 0.46, 0.51, 
and 0.56, respectively. The index of axis x, y, and z 
represent the flow, velocity gradient, and vorticity direc- 
tions, respectively. For a spherical particle, a = 4, ^ = 2. 
The temperature was determined by equilibrium calcu- 
lations before the shear flow was imposed. The temper- 
ature is ksT = 7. The initial configuration of the par- 
ticles is set to be a random distribution. A large num- 
ber of simulations were performed for volume fractions of 
0.31 < $ < 0.56 and shear rates of 5 x 10"^ < 7 < 0.1. 
These systems have particle Reynolds numbers Rep rang- 
ing from 8 X 10"'' to 1.6. 

For dispersions composed of Brownian particles in the 
steady state, the apparent shear viscosity of the disper- 
sion is defined as: 



7 



(17) 



where tr"^^ denotes the zy-components of cr"-'P'P . To ex- 
amine the system size effects in the present simulations, 
we calculated the apparent shear viscosity for three dif- 
ferent system sizes for a constant $ = 0.41, i.e., L = 32 



with Np = 50, L 



128 with Np = 3200, and L = 256 



with Np = 25600. We found only negligibly small differ- 
ences among the values of the apparent shear viscosity 
calculated for the three systems within a statistical er- 
ror. 

Figure [T] shows the dependence of the apparent shear 
viscosity on the Peclet number for several volume frac- 
tions. Here, the Peclet number is defined as Pe = 
Qurjo?'^ /ksT . For the lowest concentration ($ ~ 0.31), 
the viscosities remain nearly constant and exhibit New- 
tonian behavior. For higher concentrations ((f> > 0.41), 
the dispersions show non-Newtonian behavior. As Pe in- 
creases, the shear-thinning behavior is clearly observed 
from the higher plateau region for Pe of order 10~^ to a 
lower plateau region for Pe of about 10. From both of 
the plateau values of the viscosity curve, we can obtain 
the low shear limiting viscosity (identified as 770) and the 
high shear limiting viscosity rjoo for each volume fraction. 

In order to evaluate 770 a-nd ?7oo, we fit our simulation 
data into the following simple empirical function of Pe 
and $: 



?7/(Pe,$) = 7?o 



m - Vo 



1 + 5-i($)Pe 



(18) 



where &($) is a fitting parameter. This empirical func- 
tion is plotted in Fig. [l]with solid lines for $ > 0.46. One 
finds that it coincides with the simulation data reason- 
ably well. It is also seen that the onset of shear thinning 
appears at smaller Peclet numbers with increasing vol- 
ume fraction. We note that the particles are randomly 



distributed in the dispersions in the present simulations, 
i.e. all the viscosity data are taken in situations at which 
no shear-induced crystallization occurs. 

Figure [2] shows the dependence of the low shear limit- 
ing viscosity ?/o and the high shear limiting viscosity 7700 
on the volume fraction. Both shear limiting viscosities 
increase monotonically with the volume fraction. The 
simulation results for both of the shear limiting viscosi- 
ties agree well with the semi-empirical relations of Eq. ([2]) 
by Krieger-Dougherty, 770,00 = ?7(1 — ^/^m)"^'^*" where 
<i>m = 0.63 for 770 and <i>m = 0.73 for 7700. The dotted hue 
represents Einstein' s formula of Eq.(IT]). 

The fluid and particle inertia contribution to the ap- 
parent shear stress can be written as 



= j dxptV^Vy, 



(19) 



which is similar to the Reynolds stress of a uniform fluid 
in turbulence. This stress represents the strength of the 
hydrodynamic instability. The inertia contribution to the 



viscosity, 7y* 



'xy 



^}t/i, is about two order of 



magnitude smaller than the apparent shear viscosity. We 
conclude that the inertia contribution is negligible in the 
range 8 x 10"'^ < Rcp < 1.6. 

We next examine the dynamical motion of Brownian 
particles in a shear flow. We analyze the mean square 
displacement (MSD) in the vorticity direction (z) for the 
Brownian particles. 



1 

N, 



^(ii?i(t)-i?i(o)n. 



(20) 



Figure[n](a) shows the MSD of the dispersion with the 
lowest Peclet number for each volume fraction. These 
results are considered to be the MSD at thermal equilib- 
rium. As the volume fraction is increased, the dynamical 
behavior of the MSD varies greatly, as compared with 
that of the hydrodynamic analytical solution for a single 
Brownian particle in a shear flow [s^l • The analytical so- 
lution is derived from the generalized Langevin equation 
with memory effects. With increasing volume fraction, 
MSD increases more slowly with time. For <I> = 0.56, we 
found that a plateau region starts to appear around a 
time scale of order 10^, which is slightly greater than the 
kinematic time Tj/. This is a typical behavior of colloidal 
dispersions in glassy states, e.g., colloid glasses. 

Figure [3] (b) shows the MSD of the dispersion at the 
highest Peclet number for each volume fraction. Since 
Pe ~ 17.2 here, the shear force is much stronger than 
the thermal force. The volume fraction dependence of 
the MSD in this figure is opposite to the previous case 
at the lowest Peclet number shown in Fig. [3] (a). At 
short times, the MSD grows more rapidly in time with 
increasing volume fraction because His between particles 
are more enhanced at higher volume fraction. At later 
times, the MSDs tend to exhibit diffusive motions with 
a volume fraction independent diffusion coefhcient. 
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We then calculated the long-time diffusion coefficient 
jjsim jj-^ ^Y^e vorticity direction via 

lun^{[AR^{t)]') (21) 

t— >oo Zt 

and examined the system size effects in Df™-. In con- 
trast to the negligibly small system size effects observed 
in the apparent shear viscosity 77"^^, shows notable 

increase with increasing system size for all volume frac- 
tions. It is, however, confirmed that the diffusion coeffi- 
cients behave as 

7?f'"($, L-i) = Z?,($) - ikBT/Qin^L (22) 

with ^ ^ 2.83, where Dz{^) represents the value of diffu- 
sion coefficient extrapolated for L — > 00. The system size 
effects observed in the present simulations are essent ially 
the same as those reported in earlier papers [HI, [13, Il3l |. 
We calculated Dz{^) from extrapolation of our simula- 
tion data Dz{^,L) with finite L. Figure 4 shows the 
volume fraction and the Peclet number dependencies of 
Dz(^) normalized by = kgT /Qir-qa. For the low 
Peclet numbers, the diffusion coefficient decreases with 
increasing volume fraction, and this dependence is simi- 
lar to that of the dispersion at thermal equilibrium. On 
the other hand, for the highest Peclet number, the vol- 
ume fraction dependence of is almost constant. We 
can see that the diffusion coefficients at high volume frac- 
tions are strongly affected by shear. The dynamics of 
the Brownian particles differs considerably, depending on 
whether the thermal or shear force is dominant. The vol- 
ume fraction dependence is in qualitative agreement with 
the numerical results achieved based on the SD method, 
which is valid for long time scales [T3 |. 

Finally, the relationship between the non-Newtonian 
viscosity of the dispersions and the structural relaxation 
of the dispersed Brownian particles is examined. The 
structural relaxation time is defined as 

Tp($) = a'/Doi^) (23) 

where Dq{^) is the diffusion coefficient for each vol- 
ume fraction at zero shear. This relaxation time rep- 
resents a time needed for a particle to diffuse away 
a distance comparable to its radius. Figure [5] shows 
(jyipp _ r]^)/{r]o — T]ao) versus Tp7 for $ > 0.46. Here, 
we used -Dz($) at the lowest Peclet number as Do($), 
and rp7 can be understood as a reduced Peclet number 
Pe($) = ■ya'^ /Do{^)- One can see that the data lies on a 
single master curve 1/(1 + Ax) fairly well, where x = Tp'j 
and ^ is a fitting parameter. The shear-thinning starts 
at around Tpj ~ 1. This behavior is very similar to the 
non-Newton rheology of supercooled liquids [sH ]. 

We confirmed that the present numerical method of 
introducing the thermal fluctuation successfully repro- 
duces the fluctuation dissipation theorem for time scales 
longer than the so-called the Brownian time [2^, [2§| . Al- 
though it was confirmed also that the present method 
works quite well for the volume fractions < $ < 0.56 



considered in the present study, further careful tests must 
be needed for highly concentrated dispersions, where the 
fluctuation of a tagged particle tends to correlate with 
motions of surrounding particles. 

IV. CONCLUSIONS 

We investigated the rheological properties of monodis- 
perse concentrated dispersions with repulsive spherical 
particles by using a DNS method that accounted for a 
shear flow and thermal fluctuations. Three-dimensional 
simulations were performed at Peclet numbers ranging 
from 0.043 to 17.2 and at particle Reynolds numbers from 
8 X lO^"' to 1.6. The apparent viscosity for <I> = 0.31 
is almost constant over the Peclet number change. For 
$ > 0.41, the viscosities decrease from the plateau re- 
gion at the low Peclet number to the plateau region at 
the high Peclet number. From the viscosity versus Peclet 
number curves, we can obtain both the low and high lim- 
iting viscosities, and these results are in good agreement 
with the Krieger-Doughty relationship. The inertia con- 
tribution to the apparent shear viscosity is very small 
throughout the entire range of Peclet numbers and vol- 
ume fractions examined in the present study. As the vol- 
ume fraction increases, the behavior of the MSD in the 
vorticity direction deviates from the analytical solution 
for a single Brownian particle in a shear flow. For the 
lowest Peclet number, the MSD develops more slowly in 
time with increasing volume fraction. At <f> = 0.56, an 
onset of glassy dynamics was observed. On the other 
hand, for the highest Peclet number, the MSD develops 
more rapidly at short times with increasing volume frac- 
tion. Finally, the volume fraction dependence disappears 
at long times. The diffusion coefficient was calculated 
from the long-time behavior of the MSD in the vorticity 
direction for different Peclet numbers and volume frac- 
tions. For the lowest Peclet number, the diffusion coef- 
ficient decreases with increasing volume fraction, while 
it is almost constant over the volume fraction change 
for the highest Peclet number. We estimated the struc- 
tural relaxation time Tp from the diffusion coefficient at 
the lowest Peclet number. The present non-Newtonian 
viscosity data agrees well with a simple scaling function 
(?7°PP-r7o)/(77oo-??o) = l/(l + ATp7) for $ > 0.46, similar 
to the non-Newtonian viscosity of a model supercooled 
liquid m. 
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FIG. 1: The Peclet number Pe = Q-jTrja'^^/kBT dependence of the apparent shear viscosity of the dispersion for several volume 
fractions. The solid lines represents rff = 7700 + (?7o — '7oo)/(l + &~^(3')Pe) for $ > 0.46, where rjo is the low shear limiting 
viscosity, 7700 is the high shear limiting viscosity, and 6(4') is a fitting parameter, with b = 0.20, 0.50, and 0.63 for $ = 0.56, 
0.51, and 0.46, respectively. The dashed line represents the shear viscosity of the host fluid. 
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FIG. 2; Volume fraction dependence of the low shear limiting viscosity 770 and the high shear limiting viscosity 7700 of the 
dispersions. The solid line represents the Krieger-Dougherty relationship with $m ~ 0.63 for rjo and the dashed line <I>m ~ 0.73 
for ?7oo. The dotted line represents Einsteinfs formula of Eq.ll]). 
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FIG. 3: The mean square displacement in the vorticity direction for several volume fractions (a) at the lowest Peclet number 
(Pe = 0.0086 at $ = 0.56 and Pe = 0.022 for $ < 0.46) and (b) at the highest Peclet number (Pe = 17.2). The solid line 
represents the analytical solution for the generalized Langevin equation of a free Brownian particle [3^ . The arrows indicate 
the kinematic time — pfO^/rj — 16 and the diffusion time td = a^/Do ~ 172. 
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FIG. 4: The long-time diffusion coefficient in the vorticity direction for different Peclet numbers and volume fractions. The 
diffusion coefficient is scaled by Do = ksT /Qnria. 
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